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The Fermi energies of electrons and holes and their densities in different graphene layers (GLs) in 
the n- and p-regions of the electrically induced n-i-p junctions formed in multiple-GL structures are 
calculated both numerically and using a simplified analytical model. The reverse current associated 
with the injection of minority carriers through the n- and p-regions in the electrically-induced n-i-p 
junctions under the reverse bias is calculated as well. It is shown that in the electrically-induced n-i-p 
junctions with moderate numbers of GLs the reverse current can be substantially suppressed. Hence, 
multiple-GL structures with such n-i-p junctions can be used in different electron and optoelectron 
devices. 
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I. INTRODUCTION 



The possibility to form electrically- induced n-p and 
n-i-p junctions [l|-|3| in gated graphene layers (GLs), 
as well as lateral arrays of graphene nanoribbons and 
graphene bilayer opens up prospects of creation novel 
electronic and optoelectronic devices 0-Hj]. In contrast 
to GL structures with chemically doped n- and p-region 
in the GL structures with electrically-induced n-p and 
n-i-p junctions, there is a possibility of their voltage con- 
trol. Recent success in fabricating high-quality multiple 
GLs gEH stimulates an interest in different prospective 
devices in which the utilization of multiple-GL structures 
(stacks of disoriented GLs) instead of single-GL struc- 
tures can provide a significant improvement of their per- 
formance. Such devices include, in particular, terahertz 
tunneling transit-time oscillators (similar to that consid- 
ered in Ref . |J| ) , lasers with optical and electrical pump- 
ing, and high performance interband photodetectors pLlI — 
[l3| . Gated multiple-GL structures can be also used in 
high-frequency field-effect transistors [3] and other de- 
vices (the terahertz frequency multipliers, and plasmonic 
devices). However the penetration of the electric field 
(transverse to the GL plane) beyond the topmost GL as 
well as its sceening by electron or hole charges in GLs 
can substantially limit the influence of the gates (the ef- 
fect of the quantum capacitance (TBj). In this paper, 



we study the influence of screening in gated multiple-GL 
structures on the formation and characteristics of n- or 
p-regions and n-i-p junctions in these structures. We cal- 
culate the electron and hole Fermi energies and densities 
in GLs in the n- and p-regions as functions of the GL 
index, gate voltage, and temperature. Using these data, 
we find the voltage and temperature dependences of the 
reverse current in the n-i-p junctions with different struc- 
tural parameters. 



II. EQUATIONS OF THE MODEL 

Let us consider a multiple-GL structure with the side 
Ohmic contacts to all GLs and two split gates (iso- 
lated from GLs) on the top of this structure as shown 
in Fig. 1(a). Applying the positive (V n — V g > 0) or 
negative (V p = — V g < 0) voltage between the gate and 
the adjacent contact (gate voltage), one can obtain the 
electrically-induced n- or p-region. In the single- and 
multiple-GL structures with two split top gates under 
the voltages of different polarity, one can create lateral 
n-p or n-i-p junctions. Generally, the source-drain volt- 
age V can be applied between the side Ohmic contacts to 
GLs. Depending on the polarity of this voltage, the n-p 
and n-i-p junctions can be either direct or reverse biased. 
We assume that the potentials of the first (source) con- 
tact and the pertinent gate are ip s — and (p g — V g > 0, 
respectively, and the potentials of another gate and con- 



tact (drain) are (f g 
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FIG. 1: Schematic view of (a) a multiple-GL structure and 
its band diagrams at (b) zero bias voltage (V = 0) and (c) 
at reverse bias (for GL with k — 1). Arrows indicate the 
directions of propagation of injected electrons and holes as 
well as those thermogenerated and generated due to interband 
tunneling. 



is sufficiently wide (markedly exceeds the thickness of 
the gate layer W g separating the gates and the topmost 
GL), there are intrinsic i-regions in each GL under the 
slot. Thus the n-i-p junction is formed. The pertinent 
band diagrams are shown in Fig. 1(b) and 1(c). Since 
the side contacts are the Ohmic contacts, the electron 
Fermi energy in the k-th GL sufficiently far from the 
contacts are given by fik = ±e<^fc. Here e is the elec- 
tron charge and cpk = i p\z=kd is the potential of the fc-th 
GL, k = 1,2,..., A", where K is the number of GLs in 
the structure, d is the spacing between GLs, and the axis 
z is directed perpendicular to the GL plane with z = 
corresponding to the topmost GL and z = zk = Kd - to 
the lowest one. 

Focusing on the n-region (the p-region can be consid- 
ered in a similar way) and introducing the dimensionless 
potential tp — 2(p/V g , one can arrive at the following 
one-dimensional Poisson equation governing the poten- 
tial distribution in the z-direction (in the n-region): 



d 2 il> 
~d^ 



me 
aeVo 



K 



Y(Zl-X+).6(z-kd + d). (1) 



fc=i 



Here ae is the dielectric constant, 8(z) and S fc and 



are the equilibrium sheet densities in the fc-th GL of elec- 
trons and holes, respectively. These densities, taking into 



account the linear dispersion low for electrons and holes 
in graphene, are expressed via the electron Fermi energy 
as 



k n\hv F ) J 1 + exp (£ =F ^k/k B T) 
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where St = (tt /6)(fcsT /Hvf) 2 is the electron and hole 
density in the intrinsic graphene at the temperature T, 
vf — 10 s cm/s is the characteristic velocity of electrons 
and holes in graphene, and h and ks are the Planck and 
Boltzmann constants, respectively. Here it is assumed 
that the electron (hole) energy spectrum is e = vfP, 
where p is the absolute value of the electron momentum. 
The boundary conditions are assumed to be as follows: 
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z=-0 

Equations (1) - (3) yield 

2-Vi = r$(fa), 



z=kd+0 



0. (3) 
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for K = 1, 



W a 



for K — 2, and 
d 



jL( 2 _ fa) -fa + fa = jLr$(fa), 
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(2-fa)~fa+fa = yy-TQ(fa), 



V'fe-i - 2fa + fa +1 = — r$(V>fc), (2 < k < K - 1), 
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for K > 2. Here 
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*wo = - 



dt£ 



o 1 + exp(£ - U g tp) 
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where T = {8-7r/si)(eW g X T /V g ) cx T 2 /V g and U g 
eV g /2k B T. 
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FIG. 2: The electron Fermi energy fik as a function of the GL 
index k calculated for multiple-GL structures with different 
number of GLs K for different gate voltages V g at T — 300 K 
(upper panel) and T = 77 K (lower panel) . 



III. NUMERICAL RESULTS 

Equations (4) - (7) were solved numerically. The re- 
sults of the calculations are shown in Figs. 2-5. In these 
calculations, we assumed that se = 4, d = 0.35 nm, and 
W g = 10 nm. 

Figure 2 shows the dependences of the electron Fermi 
energy 



eV„ 



tj'h 



(8) 



as a function of the GL index k calculated for multiple- 
GL structures with different number of GLs K at differ- 
ent gate voltages and temperatures. One can see that 
the Fermi energy steeply decreases with increasing GL 
index. However, in GLs with not too large k, the Fermi 
energy is larger or about of the thermal energy. As one 
might expect, the electron Fermi energies in all GLs at 
T = 77 K are somewhat larger than at T = 300 K (see 
also Fig. 5). The obtained values of the electron Fermi en- 
ergies in topmost GLs are fii ~ 92 meV and /«i ~ 77 meV 
for V g = 1000 mV at T = 77 K and T = 300 K, respec- 
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FIG. 3: Voltage dependences of electron Fermi energies in 
some GLs in multiple-GL structure with different K at T — 
300 K. 



tively. 

Figure 3 shows the voltage dependences of the elec- 
tron Fermi energies in some GLs in multiple-GL structure 
with K = 2, K = 10, and K = 50 at T = 300 K. 

Figure 4 shows the electron densities in the struc- 
tures with different number of GLs K at different temper- 
atures. One can see that the calculated electron densities 
in GLs with sufficiently large indices (k > 15 at T — 77 K 
and T — 300 K) approach to their values in the intrinsic 
graphene (E T = 0.59 x 10 10 cm" 2 and 8.97 x 10 10 cm" 2 ). 
The electron densities in GLs in the structures with dif- 
ferent K are rather close to each other, particularly, in 
GLs with small and moderate indices. 

Figure 5 presents the Fermi energies in GLs with differ- 
ent indices at different temperatures. One can see from 
Fig. 5 (as well as from Fig. 2) that the higher T cor- 
responds to lower fik- This is due to an increasing de- 
pendence of the density of states on the energy and the 
thermal spread in the electron energies. 



IV. ANALYTICAL MODEL 

At not too low gate voltages when U g ^> 1, one can 
assume that in a number of GLs the electrons under the 
gate are degenerate, i.e., fi^ ksT, and the contribution 
of holes (nondegenerate) can be disregarded, hence, from 
Eq. (7) we obtain 



7T 



(9) 



In this case, for a single-GL structure (K — I), Eq. (4) 
yields 
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(10) 
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FIG. 4: Electron density vs GL index in multiple GL- 
structures with different number of GLs (K = 10 and K = 50) 
at different temperatures and V g = 1000 mV. 
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FIG. 5: Comparison of the vs k dependences calculated 
for different temperatures using numerical and simplified an- 
alytical (solid line) models. 



Solving Eq. (10) and considering Eq. (8), for the elec- 
tron Fermi energy in a GL (in a single-GL structure), we 
obtain 



Mi - AW 1 



(11) 



where fi g — hvp^JteVg/ 4eW g . For the same parameters 
as those used in Figs. 2- 4, ix\ ~ 150 meV. Such a value 
is markedly larger than those calculated for the topmost 
GLs in multiple-GL structures (although it is somewhat 
exaggerated because the temperature spread in the elec- 
tron energies is disregarded). This can be attributed to 
the fact that in multiple-GL structures the electron den- 
sity is shared between the topmost GL and underlying 
GLs resulting in lower Fermi energies in all of them. 

Considering multiple-GL structure with large K, one 
can neglect the discreteness of the structure and replace 



the summation in Eq. (1) by the integration. As a result, 
following Ref. [l3|, one can arrive at 



(-W g < z < 0, z > zk), (12) 



^ = lW a m (0<Z< ZK) - 
In this case, considering Eq. (7), we arrive at 

d 2 ip _ ip 2 
~dz*~l? s 

with the characteristic screening length 
L. = 



ae d 



6rc/„ 



(13) 



(14) 



(15) 



The boundary conditions for Eq. (14) are 



V>U=o=2 + W s 



dz 



dz 



0, (16) 



z=z K +0 



In multiple-GL structures with a large number of GLs 
(K ^> 1), one can extend the coordinate of the lowest 
GL to infinity and set dtp/dz\ z=oc = with ■c/'| 00 = 0. 
Solving Eq. (14) with the latter boundary conditions, we 
arrive at 



1 



(17) 



(C + z/V6L s f ' 
where C satisfies the following equation: 

C 3 - C/2 = (W g /V6L S ), (18) 

Since in reality W g >• y/6L s , one obtains C ~ 

(W g /V6L s y/ 3 oc V g 1/6 

Taking into account Eq. (8), Eq. (17) yields 



eV 



2[C+ (jfc- l)d/V6L s } 2 



Mi Ofc< 



Here 



e^s ( V g 



2/3 



(19) 



(20) 



is the Fermi energy of electrons in the topmost GL in the 
n-section (holes in the p-section), 



a k = [l + (fc-l) 7 ]- 



(21) 



and 



V6L S C W^LV 3 



V a 
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Setting d = 0.35 nm, W g — 10 nm, ae = 4, and V g = 
1 V, one can obtain L s ~ 0.44 nm. C ~ 2.14, /i! ~ 
109 meV, and 7 ~ 0.153. The [ik versus fc dependence 
obtained using our simplified analytical model, i.e., given 
by Eqs. (15) and (17) is shown by a solid line in Fig. 5. 

Equations (9) - (22) are valid when [ik > k B T, i.e., at 
sufficiently large V g or/and sufficiently small T. Since 
the Fermi energy (at fixed electron density) decreases 
with increasing temperature, the above formulas of our 
simplified (idealistic) model yield somewhat exaggerated 
values of this energy [compare the dependences in Fig. 5 
obtained numerically using Eqs. (4) - (7) and that found 
analytically using Eqs. (19)- (22)]. 




Temperature, T (K) 



V. THE REVERSE CURRENT 

The current across the n-i-p junctions under their re- 
verse bias (V < 0, see Fi g. 2 (c)) is an important charac- 
teristic of such junctions [16j . In particular, this current 
can substantially affect the performance of the terahertz 
tunneling transit-time oscillators and interband photode- 
tectors [H, [l3| • This current is associated with the ther- 
mogeneration and tunneling generation of the electron- 
hole pairs in the i-region. A significant contribution to 
this current can be provided by the injection of minor- 
ity carriers (holes in the n-region and electrons in the 
p-region). Such an injection current in the fc-th GL is 
determined by the height of the barrier for the minority 
carrier which, in turn, is determined by the Fermi energy 
(ik of the majority carrier. The latter, as shown above, 
depends on the gate voltage and the GL index. As a 
result, the reverse current can be presented as 



J = J i + J i 



Hi 



(23) 



where the injection current (which is assumed to be of 
the thermionic origin) is given by 



Ji = 



2ei>F ( k B T 



fk B TY 

\flVF J 



tt/2 



d6 cos 8 



7T /2 



FIG. 6: Temperature dependence of reverse current (injec- 
tion component, Ji) for multiple-GL structures with different 
number of GLs. 



n-region overcoming the barrier in the fc-th GL with the 
height fik is /£" ^ {1 + exp[(^ fe + VFp)/k B T}}~ 1 . Similar 
formula is valid for electrons in the p-region. The factor 
2 appears in Eq. (24) due to the contribution of both 
holes and electrons. Equation (24) is valid when the bias 
voltage is not too small: eV > k B T. At V one has 
Ji — > (as well as Jth and Jtunn)- Scattering of holes 
in the n-region and electrons in the p-region resulting in 
returning of portions of them back to the contacts leads 
to some decrease in Ji. 

The temperature dependences of the reverse current 
associated with the injection from the n- and p-region 
calculated using Eq. (24) (with the quantities [ik shown 
in Fig. 5) are presented in Fig. 6. As one might expect, 
the reverse current sharply increases with the tempera- 
ture and the number of GL. The latter is due to relatively 
low energy barriers for minority carriers in the n- and p- 
regions in GLs with large indices. For comparison, the 
injection current in a single-GL structure calculated us- 
ing Eq. (24) with Eq. (8), is given by 



J 



12J T 



d^ 



l + exp(£ + /ii/fc s T) 
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E 

k=l 



fe=i 



12J T 



K 

E 

fc=i 



1 + exp(£ + (JL k /k B T) 



1 + exp(£ + fi k /k B T) ' 



dtf 



l + cxp(e + ^/fc s T) 



(24) 



with J T = 2ev F Y, T /TT. At T = 77 - 300 K, J T ~ 0.06 - 
0.9 A/cm. Deriving Eq. (24), we have taken into account 
that the distribution function of holes which enter the 



12J T 



■ exp 



HVF\/cB Vg/SeWg 

k~B~T 



(25) 
300 K, 



At the same parameters as above and T 
Eq. (25) yields Jf ~ 0.01 A/ cm. 

When K >• 7 -1 y/fli/k^ , the main contributions to 
the reverse current is associated with GLs with large in- 
dices (in which the barriers are very low), so that one 
obtains 



Ji < KJ T . 



(26) 



As follows from Eq. (26), the injection current (and, 
therefore, the net reverse current) can be fairly large due 
to the "shortcut" by GLs with large indices (placed deep 
below the gate). 
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Since the thermogeneration is associated primarily 
with the absorption of optical phonons [17| . the perti- 
nent rate gth is independent of the electric field in the 
i-region, but it is proportional to the i-section length 21 
(I < L g ). The contribution of the thermogeneration to 
the reverse current can be presented as 



Jt 



th 



AKelgth- 



(27) 



[171 ] and < v >~ v F /2 = 5 x 10 7 cm/s [l 



113 



The quantity gth as a function of the temperature was 
calculated in Ref. [17j. Equations (24) - (27) are valid if 
21 < In, where Ir is the recombination length. In the sit- 
uations when the bias voltage between the side contacts is 
not too small (as it should be, for instance, in GL-based 
interband photodetectors), the recombination length is 
fairly long. Indeed, assuming that the recombination 
time (at T = 300 K) and the drift velocity are te 
5 x 10~ 10 

respectively, one obtains Ir ~ 250 jum. 

Assuming that 21 = 10 /an with g t h 
and g th = 10 21 cm^s" 1 at T = 77 K and T = 300 K, re- 
spectively [13], we obtain J th ~ 3.2x (lO^-lO -1 ) A/cm. 
One can see that the thermogeneration contribution to 
the reverse current is much smaller than the injection 
contribution at lower temperatures, while it can be sub- 
stantial at T = 300 K in the n-i-p structures with long 
i-region. 

The tunneling generation can significantly contribute 
to the reverse current at elevated electric fields in the i- 
region, i.e., in relatively short GL structures at elevated 
bias voltages 0, Q . This current can be calculated using 
the following formula which follows from the expression 
for the tunneling probability in GLs [l|, 0] (see, for in- 
stance @): 



evp 



eV 



ir 2 l 2 \2lhv F 



3/2 



T/3/2 

-JTJT- 



(28) 



Depending on the n-i-p junction applications, the quan- 



tities V and I should be chosen to provide either domi- 
nation of tunneling current (as in tunneling transit-time 
oscillators 0Qor its suppression (as in the interband pho- 
todetectors |13j|). 



VI. CONCLUSIONS 

We calculated the dependences of the Fermi energies 
and densities of electrons and holes in the n- and p- 
regions of the electrically induced n-i-p junctions formed 
in multiple- GL structures on the GL indices, gate volt- 
age, temperature, and the structural parameters. Using 
the obtained values of the Fermi energies and, hence, 
of the heights of potential barriers for minority carriers 
in the n- and p-regions, we found the temperature de- 
pendences of the reverse injection current for multiplc- 
GL structures with different numbers of GLs. It was 
shown that the formation of effective electrically-induced 
n- and p-regions and n-i-p junctions, i.e., the n-i-p junc- 
tions with suppressed reverse currents in multiple- GL 
structures with several GLs is possible. The utilization 
of the electrically-induced n-i-p junctions in multiple-GL 
structures in different devices, such as terahertz tunnel- 
ing transit-time oscillators, lasers, high performance in- 
terband photodetectors, and some others might provide 
an enhancement of the device performance (an increase 
in output power and responsivity) and widening of their 
functionality (owing to the possibility of the gate- voltage 
control) . 
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